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We observe the formation of heaps of repulsive microspheres, created by flowing a colloidal mi- 
crosphere suspension towards a flat-topped ridge placed within a quasi two-dimensional microfluidic 
channel. This configuration allows for both shear and normal forces on the microspheres in contact 
with the ridge. The heaps, which form against the ridge, are characterized by two distinct phases: 
a solid-like bulk phase in the interior, and a highly- fluctuating, liquid-like state which exists along 
its leading edge. We observe that heaps only form above a critical flow velocity, v c , and that they 
are destroyed by thermal rearrangements when the flow ceases. We monitor the dynamics of heap 
formation using fluorescence video microscopy, measuring the heap volume and the angle of repose 
in response to microsphere deposition and erosion processes. We find that the steady-state angle of 
repose, Of, increases as a function of inflow velocity, Voo, with a functional form Of oc V^oo — v c . 



I. INTRODUCTION 

Flowing particles will accumulate against an obstacle 
in their path, a phenomenon exhibited by industrial or 
agricultural grains accumulating in conical piles and mi- 
croscopic particles undergoing membrane filtration, par- 
ticularly ultrafiltration*^ This similarity in behaviors has 
raised the question of whether there is an underlying 
commonality which explains the general phenomenon, 
despite the large difference in length scales, the vastly 
different importance of Brownian motion and kinetic en- 
ergy, and the near-absence of interparticle friction in col- 
loidal interactions. 

The quest to understand the formation of piles of gran- 
ular materials has, in fact, become a model system in 
which to investigate how bulk behaviors arise from highly 
cooperative effects at the grain scale. Although the ba- 
sic mechanism for understanding the angle of repose has 
long been knownj^ interest in sandpiles was rekindled by 
the pioneering work of Bak et alM. Since then, detailed 
studies of pile-formation have led to a quantitative un- 
derstanding of the intermittency of granular dynamics 
near the threshold of flowp^ force chains, 6 7 how local 
dynamic rearrangements^"^ contribute to pile instabil- 
ity, the prediction of avalanches mem ory effects,^ fric- 
tional rheology in submerged flows pHUl and the relation- 
ship between internal friction and the cessation of flowJ^ 

In this paper, we report on the formation of sandpile- 
like structures in a two-dimensional (2D), approximately 
hard-sphere, colloidal system by flowing a dilute suspen- 
sion against a ridge within a microfluidic channel. Un- 
like granular piles that are formed by gravity, the flowing 
colloidal suspension is driven by a static pressure gradi- 
ent which partitions the system into coexisting gas-like, 
liquid-like, and solid-like microsphere phases. We find 
that for sufficiently large Peclet number (corresponding 
to large flow velocities), stable heaps of microspheres are 
able to form by the deposition of particles out of a di- 
lute suspension. When the flow is reduced, the heap 
re-evaporates in a thermally activated, reversible pro- 
cess. Below, we quantify the conditions under which 
these heaps form, and report their key properties as a 



function of the flow velocity. 

Our work extends earlier work on dense 2D collo idal 
systems in which a complex phase diagram ] 16 * 17 ! the 
influence of hard and soft interparticle potentials* 18 * 19 * 
the onset of layering*^*, the dynamic response to 
supercooling*^ and the effect of confinement!^* were 
demonstrated. Analogous granular experiments showed 
a phase space that included even a colloid-like ordered, 
fluctuating phaseP^Here we add the action of a combined 
compressive and shear field to the classical 2D colloid. In 
dense bulk suspensions the application of similar com- 
pressive gravitational fields^* 2 ^* anc [ shearing fiel d 26 * 27 * 
have driven gas-fluid, fluid-solid and solid-solid transi- 
tions. Our system joins a growing range of techniques 
that utilize the tunability of key suspension parameters^* 
(Peclet number, polydispersity, interaction potential) to 
explore the interplay between Brownian, steric, and spe- 
cific interparticle potential effects. 



II. EXPERIMENTAL SETUP 

We perform experiments with a dilute suspension of 
fluorescent d = 520 nm (5.9% polydispersity) polystyrene 
microspheres (Bangs Laboratories) at concentration 
0.1 % (w/w), corresponding to 6.7 x 10 11 micro- 
spheres/mL, suspended in a buffer solution described 
below. The suspension is pumped through a shallow 
microfluidic device with a flat-topped ridge rising from 
the channel floor. The full length of the microchannel is 
14 mm, long enough to hydrodynamically decouple the 
access holes from the ridge. The ridge runs perpendicular 
to the mean flow direction so that microspheres accumu- 
late upstream of it (see Fig. [I]). The device channel has 
a quasi-2D geometry, with width W = 200 d and height 
H = 1.8 d. The ridge has width w = 100 d and occu- 
pies the central 50 % of the channel width, permitting 
flow around the sides, and a height h = 1.5 permitting 
flow over the top. Because d > H — h, the streamlines 
which intersect the ridge push microspheres against it 
with both normal and shear forces(see Fig. [l]b and c). 
Crucially, the geometry is resilient to catastrophic clog- 
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FIG. 1. (a) Schematic of microchip geometry (not to scale), 
showing a single microsphere flowing towards the flat-topped 
ridge over which fluid can flow. (b,c) Raw fluorescence images 
with linear grayscale (black = background, white = highest 
fluorescence intensity) of heap at steady state, at areal density 
p — 235/(100/im) 2 and inflow velocity Voo — 33 pm/s. In the 
two seconds between (b) and (c), three tracked microspheres 
marked by the circle, square, and triangle symbols are shown 
in both figures to illustrate their relative displacements. 



ging since neither the ridge nor the resulting microsphere 
heaps span the entire width micro channel. 

The microfluidic channels are manufactured on a sil- 
icon wafer, using photolithography to define the chan- 
nels and reactive ion etching to transfer the pattern into 
the wafer. Access holes « 300 /im in diameter are sand- 
blasted through the silicon wafer with 17 pm alumina 
particles. The fluid system is closed by bonding the sil- 
icon wafer to a borosilicate coverglass (Corning) coated 
with a 20 pm layer of cross-linked PDMS. This hybrid 
PDMS-glass coverglass allows for reversible chip bond- 
ing, allowing many experiments to be performed on the 
same chipP^ESl The finished device is mounted to a sam- 
ple holder that provides 10-/iL reservoirs to supply sus- 
pensions to the access holes. The suspension flow rate 
through the channel is set by a digital pressure regu- 
lator (AirCom PRE1-UA1) that applies pressurized air 
(0 — 10 kPa) into the o-ring sealed reservoir. 

To introduce the suspension into the chip without mi- 
crospheres interacting with air- water interfaces, we first 
wet the chip with a microsphere- free buffer. Addition of 
1 pM sulforhodamine dye to the buffer allows verification 
of complete wetting of the devices by fluorescence mi- 
croscopy, especially of the gap above the ridge. We then 
remove ~ 90% of the wetting buffer from the reservoir 
and inject the microsphere-bearing suspension into the 
remaining buffer. This procedure ensures air bubbles are 
not introduced into the access hole, and thereby allows 
us to maintain a stable flow rate over the duration of the 
experiments. At the start of each experimental run, we 
flow the suspension into the device. 

Using fluorescence microscopy, we both record the lo- 



cal microsphere distribution at the ridge, and continu- 
ously determine the two free suspension parameters: in- 
flow velocity and mean 2D areal density p, given 
by the number of microspheres per unit area. To do 
this, we use microspheres that are fluorescently labeled 
with DragonGreen dye (480 nm excitation, 520 nm emis- 
sion, Bangs Laboratories). We obtain the inflow veloc- 
ity Voo by finding the angle of maximum autocorrelation 
in two-dimensional graphs of intensity as a function of 
time and vertical position for each image column. In the 
Hele-Shaw geometry, the lateral velocity field far from 
the heap is spatially uniform,^ except within a distance 
H of the channel edges. Thus measurement of over 
47 pmx 100 /im centered 200 /im upstream from the ridge 
suffices to fully characterize the bulk flow. We use the 
same spatial region to measure p by counting the number 
of extended maxima above a noise threshold. 

Experiments are performed on a Nikon Eclipse 80i fluo- 
rescence microscope, where the microsphere fluorescence 
is excited by a continuous- wave broadband Xenon lamp 
(X-cite 120). Emitted light is collected with a 10X Plan 
Fluor objective, and 658 x 496 /im 2 images are acquired 
with an Andor Luca EMCCD camera. All images pre- 
sented are oriented such that the mean inflow and out- 
flow directions are vertical. The physical pixel size is 
10 x 10 /im 2 , while the diffraction limited projection of a 
point object on the CCD at our emission wavelength is 
8x8 /im 2 . This is larger than an equivalent microsphere 
diameter, and we are thus not able to resolve individual 
microspheres in dense piles. The total depth of field is 
5.6 /im, larger than the depth of the microchannel. In the 
absence of external driving, we measure the diffusion con- 
stant of a dilute set of microspheres in the microchannel, 
far from the ridge, to be D = 0.55 pm 2 / s. Therefore, the 
microsphere diffusive time scale is td = d 2 /D = 0.49 s in 
the microchannel, so we set our camera's exposure time 
to r e — 0.1s. This exposure time is sufficiently short 
to observe microsphere thermal rearrangements because 
T e < Tp. As we will discuss below, we measured a lowered 
and non-constant diffusion coefficient due to wall-dra^^ 
which lead to increased actual m- 

For the microspheres to have reversible contacts, it 
is necessary to reduce attractive interactions. We use 
a combination of steric and electrostatic stabilization. 
For steric stabilization, we add a non-ionic surfactant 
0.1 % (v/v), equivalently 1.6 mM, Triton X-100 (Sigma 
Aldrich) to a stock suspension of concentration 0.1 %w, 
which is enough surfactant to reach full coverage of the 
microspheres if the surface area per surfactant molecule 
is approximately 12 A 2 . To ensure that the surfactant 
coverage of the microsphere surface approaches full cov- 
erage, we tumble the suspension over 2 hours at « 10 rpm 
to disrupt the surfactant micelles (CMC for Triton X-100 
w 0.02 % (v/v), equivalently 0.22-0.24 mM). To test that 
we have successfully disrupted the micelles into covering 
the microspheres, we measure the particle size distribu- 
tion of the tumbled solution using dynamic light scatter- 
ing. While the presence of free surfactant introduces an 
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entropic attractive force, the depletion force, we observe 
that its magnitude is not sufficient to form weakly aggre- 
gated gels over months in the absence of fluid flow, much 
less in the presence of flow. Thus, for our purposes, this 
depletion force is negligible. 

We introduce a long-range repulsive electrostatic inter- 
particle interaction by using microspheres with sulfate- 
modified surfaces (pKa ~ 2.5). The interaction strength 
can be tuned via the buffer component of the suspension. 
We set pH = 5.4 with a 10 mM citric acid buffer in water 
titrated with sodium hydroxide. Using the Henderson- 
Hasselbalch equation, we estimate the concentration of 
free ions and calculate a Debye length of 3nm. Exper- 
imentally, the ionic strength of the solution is adjusted 
by dilution of a high-concentration buffer with deionized 
water until the conductivity is 0.2 ±0.02 mS/cm. The mi- 
crosphere ("-potential is ~ —43 mV, ensuring the electro- 
static stabilization of the suspension.^ All conductivity, 
^-potential, and DLS measurements are measured with a 
Zetasizer (Malvern). 



III. RESULTS 

An illustrative image of a heap at steady state is shown 
in Fig. [TJdc. The inflow (top to bottom) carries micro- 
spheres towards the ridge that runs horizontally. Above 
the ridge, a dense, approximately triangular zone forms. 
We measure the geometry of the dense region by finding 
the points of steepest intensity gradient for each pixel 
column of the image, and taking the distance between 
them as the local extent e(x,t) for each column, where 
x is the horizontal position and t is the time. From each 
e(x,i) profile, we compute the heap area by integrating 
A(t) = J e(x,t)dx. An angle of repose 6{t) was defined 
from the slope of e(x, t) at the outer edge of the ridge. 

Although the suspension flow pattern around the ridge 
is disturbed by its presence, the microsphere-excluded 
zone below the ridge provides a clear indication of 
amount of liquid passing through the dense zone and over 
the barrier. We characterize the magnitude of this flow 
by the asymptotic physical extent of the excluded zone 
£ (see Fig. [T]o) . In Fig. [l]3c, the high microsphere areal 
density indicated by the bright bands delimiting the ex- 
cluded zone is primarily the result of deterministic lateral 
displacement!^ of the streamlines colliding with the heap. 
While the rate of microsphere erosion from the heap also 
contributes to an enhanced p, this rate is small at steady 
state. 

We varied the incoming flow velocity over a range of 
^oo = 1 — 100 fim/s to explore the nature of the dense 
heap formed against the ridge. In Fig. |2j we plot the 
pixel-level fluctuation amplitude of the fluorescence in- 
tensity over 0.8 s (80 frames) after a steady-state heap 
was formed (see below for definition of steady state). Un- 
der this analysis, microsphere-free regions and static mi- 
crospheres appear dark. Moving, isolated microspheres 
appear as streaks which blend to a homogeneous gray at 




FIG. 2. Maps of mean density fluctuation amplitude of mi- 
crosphere heaps formed at inflow velocities (a) Voo = 2 /xm/s, 
(b) Voo — 6 jum/s, (c) Voo — 22 /mi/s, and (d) Voo — 84 /xm/s. 
The bulk flow is from top to bottom. Each pixel is given a 
grayscale value proportional to the standard deviation mea- 
sured over 8 s (80 images) after steady state was reached. 
Images in the left column have the equal horizontal and ver- 
tical scales, images in the right column are stretched by a 
factor of 10 in the vertical direction. Bright areas at the sur- 
face of each heap are strongly fluctuating (fluid- like), while 
dark areas have a low fluctuation amplitude (solid- like). 
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FIG. 3. Time-series of (a) microsphere concentration p, (b) 
inflow velocity ^oo and scaled width of excluded zone £,/W (c) 
heap area A, and (d) angle of repose in response to a six 
hour pulse of increased velocity. See Electronic Supplemen- 
tary Information for a video of these dynamics. 



high flow velocities. The dense region appears gray, and 
the gray scale indicates the local magnitude of micro- 
sphere rearrangements (without correction for the local 
p). In all images, we observe a highly fluctuating bright 
layer on the upstream side of the dense particle region. 
For sufficiently high flow velocities, we observe a region 
of low fluctuation amplitude between the ridge and the 
bright layer (Fig. [2bcd). This region is not present at 
low flow velocities (Fig. |2^i) . We interpret the presence of 
the low-fluctuating region as the formation of a quasi-2D 
heap, by analogy with granular heaps. In the following, 
we will refer to the region of low fluctuation amplitude 
as solid-like and the region of high fluctuation amplitude 
as fluid-like. Note that the fluid-like phase appears to be 
compressed by the applied flow field in a similar way that 
an atmosphere can be compressed, but that the large p 
makes it unlikely that the particles can be treated as a 
classical gas. 

We observe that the heaps in Fig. [2] do not form per- 
fect triangles, even in the steady state. Instead, we find 
heaps with rounded peaks and convex sides at low 
(Fig. |2}dc), and pointed tips and concave sides at high 
(Fig. |2]i). We also observed a small degree of asym- 
metry between the left and right side of the heap. 



To define the process and time scale over which the 
heaps reach steady-state, we perform experiments in 
which we form heaps with solid-like cores at constant , 
and subsequently drop Voo to a value at which only the 
fluid-like layer was previously observed. The results are 
shown in Fig. [3] Initially, the device was primed with 
suspension at = 2/im/s for 3 hours without heap 
formation. After Voo was quickly increased to 84/xm/s, 
we observed an asymptotic increase of the heap area to 
A « 6000 pm 2 and a steady-state angle of repose of 
Of = 6.2°. The final heap contained ~ 2 x 10 4 parti- 
cles and was ~ 56 <i tall at its peak. When was re- 
turned to 2 jitm/s, a short-lived expansion of the heap was 
recorded. Thereafter we observed a gradual reduction of 
A and until only the fluid-like layer remained. The 
process thus appears as if the solid-like phase gradually 
erodes, in a process akin to evaporation or melting. All 
stages of the experiment indicated that Brownian motion 
of microspheres is essential: an atmosphere is present on 
the forming heap at all times, the heap impulsively ex- 
pands when is reduced, and the heap gradually re- 
duces thereafter. Visual inspection of the steady state 
dynamics and of the heap formation and erosion process 
are available through the Electronic Supplementary In- 
formation. 

For this experimental run, we determined a deposition 
time constant of 42 min and an erosion time constant of 
17 min by fitting exponential relaxation laws to A(t) at 
early times (Fig. [3^); similar values are obtained for 0(t). 
At longer times, we observe a non-exponential, grad- 
ual increase of A(t) and 0(t); we believe this may result 
from the slow increase of the microsphere areal density 
throughout the high-v^ phase from p = 50/(100/im) 2 to 
nearly twice that. Note that this corresponds to a change 
in the interparticle spacing of the suspension from 16 d 
to lid, indicating that we remain in a weakly interact- 
ing gas-like flow regime during this time. We observed 
that variations within this range did not have a strong 
effect on the steady state, but did not conduct a detailed 
investigation of the effect of p on the dynamics or final 
state. 

Fig.[4]shows a detailed exploration of the functional re- 
lationship between heap geometry and , as illustrated 
by the images shown in Fig. [2j At each value of v^, 
we allowed the heap to reach a steady state, defined by 
a change in A and smaller than 1%/hour (durations 
ranged from 2 — 6 hours). To distinguish the cross-over 
to a heap- forming state, we observe that the two exper- 
iments with Of « 0.5° are unambiguously heap-forming 
based on visual identification of the characteristic low- 
fluctuation amplitude region using the temporal standard 
deviation technique. In this regime, is difficult to mea- 
sure accurately due to the small height of the heap (~ 2 
pixels). 

We observe that there is a well-defined relationship be- 
tween and Of, with finite-sized piles forming for values 
of above a critical velocity v c ~ 5/im/s. Note that 
Fig. [4] is derived from two distinct experimental series 
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FIG. 4. Steady state angle of repose ^/ as a function of in- 
flow velocity Voo- Solid points denote experiments in which a 
heap was formed and an angle could be measured; open sym- 
bols denote experiments in which no heap (NH) was formed. 
Velocity error bars represent our measurement uncertainty 
in the mean measured, while the angle error bars represent 
the min/max range of angles observed at steady state. Red 
points correspond to the increase/decrease sweep shown in 
Fig. [5j with A indicating the runs in which Voo was increasing 
and T indicating decreasing. Inset: All solid (heap- forming) 
points plotted on logarithmic axes as a function of Voo — v c , 
where v c = 5/xm/s. Dashed line provides a comparison to a 
power law with exponent 1/2. 



with two different instances of device assembly, and that 
p varied but was kept within the dilute regime. For the 
black points in Fig. [ij p ranges from 40 to 150 //im 2 , 
while for the red points, p ranges from 290 to 370//im 2 . 
The relative order of values was changed between the 
runs, and the possibility of hysteresis was checked via a 
series of pressure changes in which was systematically 
increased and then decreased to the same values. This 
case will be discussed in more detail below. We observe 
that the data, in spite of the variation in p, approximately 
collapse to a single curve which take the approximate 
form Of oc \/voo — v c . 

The detailed results of the hysteresis-probing experi- 
ment are shown in Fig. |5j The applied pressure is in- 
creased and decreased along a staircase profile, resulting 
in a velocity step of 6 pm/s every two hours (see Fig.[5^i), 
sufficient to reach steady state. The resulting A(t) and 
6(t) are shown in Fig. [5]3c, and the triangular points in 
Fig. [4] are derived from these plots. We observe only 
weak hysteresis, indicating that the steady state geom- 
etry is independent of whether driving is increasing or 
decreasing. We further note that the apparent hysteresis 
is correlated to the variation in p. 

An interesting feature of Fig. [5] is the presence of sud- 
den variations in 6 (see Fig.^), which are not present in 
the pile area curve (Fig. [5^). Visual inspection of the raw 
video reveals that these likely represent rearrangement 
events on the heap surface that affect only the profile 



l zl 

8 
> 



30 
20 
10 



(a) 




FIG. 5. Time-series of (a) inflow velocity Voo, (b) heap area 
A, (c) angle of repose for a series of increasing and de- 
creasing pressure steps. The areal density for this run is 
p = (330 =b 42)/(100/xm) 2 , where the range covers two stan- 
dard deviations. 



shape, but not the profile area. Interestingly, these rear- 
rangements are more prominent during the second half 
of the experiment (erosion) than during the first half (de- 
position). 



IV. DISCUSSION 

There are significant similarities between granular 
heaps created by the raining of particles under grav- 
ity, and the microsphere heaps created in the present 
study: both are largely triangular in cross-section and 
can be built with different angles of repose. However, 
granular materials have an angle of repose which is de- 
termined by the static friction coefficient p of the parti- 
cle material,^ with larger angles supporting greater tan- 
gential forces through larger values of p. For example, 
Lemieux and Durian^ observed « 20° (equivalent to 
p — tan# = 0.36) for heaps formed from dry, millimeter- 
scale, spherical, glass beads. 

Compared to granular piles, we measure much lower 
angles of repose from 0.5° to 7° (see Fig. |4|, and have 
no expectation of frictional contacts between the micro- 
spheres. Therefore, a closer comparison lies with simula- 
tions and experiments on frictionless particles. Peyneau 
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and Roux 40 observed lower limit of ~ 5.8° for a molec- 
ular dynamics simulations of a frictionless pack under 
shear, independent of the applied pressure. Similarly, 
for 120 jam bubbles rising against to form an upside- 
down heap, Lespiat et al. 41 observed a lower limit of 
« 4.6°. While these values are consistent with our 
observations, and support the expectation that friction 
play a minimal role, we find that we are able to form 
steady-state heaps with even lower angles of repose, ap- 
proaching Of « 0° (see Fig. [3]). This suggests that at 
low Peclet number, Brownian microspheres are able to 
explore their cages sufficiently wel l 33 ! 43 ! to destroy the 
geometrical constraints which give rise to the finite angle 
of repose in frictionless granular/bubble systems. The 
fluctuations also lead to lubrication and frictionless con- 
tacts since even small thermal displacements are able to 
disrupt the contact network.^1 

We have described the microsphere heaps formed 
above v c as having a solid-like core, which we identified by 
its small temporal density fluctuations when compared to 
the dense region formed below v c , or the boundary layer 
between the heap and the surrounding dilute suspension 
(Fig. [2]). Whether the microspheres exist in a solid-like 
or liquid-like state depends on a competition between 
the stabilizing normal component of the drag force of the 
fluid flowing through the heap, its de-stabilizing tangen- 
tial component, and the de-stabilizing Brownian motion 
of microspheres. At the lowest inflow velocities (such as 
in Fig. [2^1), the Brownian motion dominates and we ob- 
serve only a liquid-like phase, without forming a heap. 
The liquid-like region of our system mirrors the density 
increase (concentration polarization) that is observed at 
the interface of a filter membrane in cross-flow filtration 
systems.^ At sufficiently high suspension concentrations 
or flow velocity through the membrane, these systems 
also show the formation of a filter cake, mirroring the 
solid-like bulk phase of the interior of the heap that we 
observe at high (Fig. J2j)cd). Note that cross-flow fil- 
tration differs from our experiment in that the ratio of 
normal to tangential flow components is far larger in our 
system. It is commonly assumed that the cake phase is 
solid per se, while our heap still undergoes thermal reor- 
ganization. 

Song and Elimelech 45 describe a similar fluid-like layer 
as having an exponential density distribution. They de- 
fined a criterion for cake formation based on the filtra- 
tion number Np, defined as the ratio of the energy re- 
quired to bring a particle from the surface of the cake 
to the interface of bulk solution and concentration po- 
larized layer. Our measurements support a similar cri- 
terion, in which we examine the relative importance of 
Brownian motion and fluid drag for the transition from 
a solely-liquid-like state (Fig. [2k) to one in which both 
phases are present (Fig. J^cd). We quantify the im- 
portance of Brownian motion with the Peclet number 
Pe = to/to,, where r a = d/v is the particle advection 
time and tjj = d 2 / D is the self-diffusion constant D. We 
measured D = 0.25 /am 2 / s in a dense suspension that was 



created by ceasing the flow around a liquid-like state. A 
simple argument demanding flow continuity through the 
heap and over the ridge yields that the mean normal 
velocity of the aqueous suspension at the shear bound- 
ary for particles is ^heap — v oo^- Using our measured 
value of D, we obtain 1 < Pe < 23 in the incoming di- 
lute suspension flow for the range of explored here. 
The critical velocity v c ~ 5 /im/s corresponds to Pe « 2. 
This argument is only an estimate, as it neglects the influ- 
ence of both shear and the incident particle density, both 
of which likely influence Of. In addition, the diffusion 
constant is modified by geometry- and concentration- 
dependent contributions to the S tokes drag j 46 | 4T | as well 
as the degree of confinement ) 32 l 43 l Nonetheless, in exper- 
iments on dense suspension flows, Isa et al.^i also ob- 
served a change in rheological properties at similar values 
ofPe. 

The microfluidically- assembled microsphere heaps will 
also provide a convenient system in which to examine 
the behavior of a dense, solid-like colloidal system as a 
function of the applied stress. In Figs. [3] and [5] we note 
a transient compression and expansion of the heap upon 
sudden increase and decrease of the compressive stress 
set by ^oo, respectively. For instance, at t = 533 min 
m Fig. the area heap transiently expands by 18% 
under a drop in pressure of « 30 Pa, which corresponds 
to a modulus of ~ 10 2 Pa, far smaller than the few 10 9 
Pa typical for polystyrene. This observation provides ev- 
idence that individual microspheres are not in constant 
contact with each other but instead continually exploring 
cages by Brownian motion. 

Finally we note that the experimental system pre- 
sented here lends itself for repeated analysis of non- 
steady state systems that may be driven into super- 
critical regimes. For instance, we are able to create 
patches of colloids close to space-filling packing by re- 
leasing close-packed heaps through sudden stopping of 
the flow. As such, this system will be useful for explor- 
ing the temperature-dependence of jammed solids near 
the jamming transition! 49 * 50 * Our observation of a crit- 
ical v c above which jamming occurs is consistent with 
suggestions^ that thermal systems should jam by in- 
creasing the osmotic or hydrostatic pressure. While not 
explored in the present study, the interparticle potential 
is tunable by adjusting the pH of the buffer component 
of the suspension. 



V. CONCLUSION 

In this Paper, we have established a new system for 
forming a solid-like heap of microspheres within a mi- 
crofluidic flow. We observe that the inflow velocity is 
the key control parameter in setting whether a steady- 
state heap is formed (large v^, small Peclet number) 
or destroyed (small v^, large Peclet number). The 
steady-state angle of repose of this heap increases from 
Of « 0° near the transition, to larger angles (up to 
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7°) at the fastest velocities explored, with a functional 
form Of oc v^oo — v c . We propose that this effect arises 
through the suppression of thermal fluctuations at larger 
flow velocities through reaching of a critical density of 
microsphere in a region of concentration polarization. 
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